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Abstract. We consider a two-dimensional massless Dirac operator H in the 
presence of a perturbed homogeneous magnetic field B = Bo + b and a scalar 
electric potential V. For V £ Lf oc (]R 2 ), p £ (2, do], and b e L^ oc (IR 2 ), q £ 
(1, oo], both decaying at infinity, we show that states in the discrete spectrum 
of H are superexponentially localized. We establish the existence of such 
states between the zeroth and the first Landau level assuming that V = 0. 
In addition, under the condition that b is rotationally symmetric and that V 
satisfies certain analyticity condition on the angular variable, we show that 
states belonging to the discrete spectrum of H are Gaussian-like localized. 

1. Introduction 

Graphene is a two-dimensional lattice of carbon atoms arranged on a honeycomb 
structure. Due to its unusual properties it has attracted a great deal of attention 
since its discovery [4] [2"T| . One of the striking facts about graphene is that the 
dynamics of its low-energy excitations (the charge carriers) can be described by 
massless two-dimensional Dirac operators. An interesting feature of Dirac fermions 
is the lack of localization under the influence of an external electric potential [30jll5j. 
This fact, related to Klein's paradox [4], is due to the peculiar cone-like gapless 
structure of the spectrum of massless free Dirac operators. 

It was suggested in [7] that it is possible to confine such Dirac fermions in 
graphene by inhomogeneous magnetic fields of the type B = Bq + b, where Bo > 
is a constant and b a perturbation with negative flux that decays at infinity. The 
spectrum of the corresponding Dirac operator in a constant magnetic field Bq is 
given by the (relativistic) Landau levels. The idea is that as the perturbation b 
is turned on eigenvalues will emerge from the Landau levels giving rise to states 
localized on the bulk of the support of b. In this manner a so-called (magnetic) 
quantum dot or artificial atom can be created. These type of models, also with an 
external electric potential V, have been further studied in the physics literature, for 
instance, in [5J [331 133 ES] f° r the one particle case and in [131 H] for the multiparticle 
case. The articles 8, 22] [3TJ [TB] deal with specific electromagnetic fields for which 
the model is partly solvable or suitable for numerical computations. 

In this article we consider a large class of electromagnetic perturbations (b, V) 
with V G Lf oc (]R 2 ), p G (2,oo], and b G L? oc (K 2 ), q G (l,oo], both decaying at 
infinity. The essential spectrum of the corresponding Dirac-operator H describing 
the quantum dot is given by the Landau levels. We show that eigenfunctions 
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belonging to the discrete spectrum of H are superexponentially localized, i.e., they 
decay faster than any exponential. In the case when V = we verify the existence 
of eigenvalues between the zeroth and the first Landau-level assuming that b < 
0. Assuming that a certain analyticity conditions on the angular variable of V 
is fulfilled and that b is rotationally symmetric we prove that those states are 
actually Gaussian-like localized. These type of results on superexponential and 
Gaussian localization, although new for Dirac operators, are known to hold for 
spinless magnetic Schrodinger operators [SJ [TUl HH1 HE] • We benefit from this insight 
to prove our statements. A precise description of our results is given in the next 
section. 



2. Results 

We consider the massless two-dimensional Dirac operator with an external mag- 
netic field B : R 2 —> R, pointing perpendicularly to the plane, and an electric 
potential V : R 2 -> K. We are interested in the Hamiltonians 

(1) D A := a ■ (p - A), 

(2) H:=D A + V, 

a priori defined on C£°(R 2 ;C 2 ) C L 2 (R 2 ;C 2 ). Here p := fV is the momentum of 
the particle and <x := (01,(72) is a vector whose entries 

1 \ / -i 



,T| 1 i o )■ CT2 - v i o 

are Pauli matrices. The magnetic field B enters in the definitions ([IJ and © by 
means of the magnetic vector potential A = (Ai,^) : R 2 — > R 2 through the 
relation 

(3) B = d±A 2 - d 2 A 1 =: curl A , 

which is understood in the sense of distributions. 

Throughout this article we assume the following on (B, V): 

(Al) B = Bo + b where B > is a number and b e L|' oc (R 2 ;R) for some 
q € (l,oo] and lim^oo || l{| x |> n }6||oo = 0. 

(A2) V € L p loc (R 2 ;R) for some p £ (2,oo] and linwoo ||l{|x|>n}V||oo = 0. 
Here !/(•) denotes the characteristic function on the set /. Assuming that B fulfills 
(Al) we can always find A 6 L[ oc (R 2 ;R 2 ) for some t £ (2,oo] satisfying ([3]) (see 
Remark [7]) ■ For such magnetic vector potentials and electric potentials V satisfying 
(A2) we know that the operators defined in (fTJ) and © are essentially self-adjoint 
(see Subsection 13. ip . We denote their self-adjoint extensions by the same symbols 
and their domains by X>(Da) and T>(H) respectively. 

To the homogeneous magnetic field Bq we associate the vector potential 

(4) A := -j-(-x 2 ,x 1 ), 

satisfying curlAo = Bq. It is well known that the spectrum of -Da consists of 
infinitely degenerated eigenvalues (l n )nez, called Landau levels, given by 



l n ■= sgn(n) y/2\n\B , n £ Z , 
where sgn(n) = n/\n\ if n ^ and equals one if n = 0. 
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Given a self-adjoint operator T we write cr pp (T) , Cd (T) , and a CS s(T) to denote 
the pure point, discrete, and essential spectra of T respectively. Our first main 
result is as follows. 

Theorem 1. Assume that B satisfies (Al) and let A € L\ Qc (R 2 ; M 2 ) , p € (2, oo], 
with curl A = B. Then, the spectrum of D A is symmetric with respect to zero and 

Cess(-DA) = (ln)neZ ■ 

Moreover, 

(a) If b < and strictly negative on some open set, then the discrete spectrum 
of Da on (0, Zi) is non-empty, i.e., ad(D A ) H (0, Zi) ^ and 

dim(Raa(l( ,/ 1 )(I>A))) = oo . 

(b) Ifb > then 

dim(Ran(l ( o, il )(£> A ))) =0. 

This theorem is a consequence of lemmas [2] and [3J That the spectrum of D A is 
symmetric with respect to zero is well known; see, however, Proposition [T] 

Remark 1. A similar result to Theorem\7\is shown in [3] when b is replaced by 
Xb and A is assumed to be sufficiently large. Moreover, in [3] stronger regularity 
assumptions on b are made. In addition, the magnetic vector potential a associated 
to b is assumed to decay at infinity. However, the results of [3] hold for more general 
background magnetic fields than Bq . We note also that our proof differs from the 
one in [3] . 

Remark 2. Assume that (Al) and (A2) are fulfilled. As a conseguence of lemmas 
[7] and\^ below, 

a ess (H) — a css (D A ) — a css (D Ao ) , 

for any A € Lf oc (R 2 ;K 2 ), p € (2,oo], with curl A = B. 

Our next result state that eigenfunctions corresponding to the discrete spectrum 
of H are super-exponentially localized. 

Theorem 2. Assume that B and V satisfy (Al) and (A2) respectively and let A £ 
iJ 5 QC (M 2 ; M 2 ), p € (2, oo], with curl A = B = Bq + b. Then, for any eigenfunction ^ 
of H = D A + V with iJ* = and E € M \ o-(D Ao ) the following holds: 
For every r € [2, oo] and 7 > there exists an R > such that 

(5) ||l { | x |> fl . }e ^ x lvl/|| r <oo. 

This theorem is proven in Section [5] 

Remark 3. This type of results are known to hold for magnetic Schrddinger op- 
erators (p — A) 2 + B. Our proof follows the ideas presented in In fact, since 
our operator is linear in A, some parts of the argument are more straightforward. 
For instance, we do not require that b € C 1 (M 2 ;M) decays in the C x -norm as done 
in [B]. 

Remark 4. One essential ingredient in the proof of Theorem fj| is the explicit 
knowledge of the Green function Go of D Aq . This is calculated in AvvendixYfi\ 

In order to obtain Gaussian decay we make further assumptions on (B, V). Let 
T = M/(2ttZ) and let v = v{r,9), (r,9) € M+ x T be the potential V written in 
polar coordinates. We assume: 
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(A3) B is radially symmetric, i.e., &(x) = b(r), r = |x|. 

(A4) For any (r, 6) e M+ x T the mapping K 3 a h-> u(r, + a) =: v Q (r, 0) has an 
analytic continuation v z (r,9) to C. Moreover, for any r > there exist a 
p £ (2, oo] and a real- valued function u r € Lj ) oc (R + x T,rdrd6) such that 
II l{r>n} u r||oo — ^ as 7i — y oo and 

|«,(r,0)| <tir(r,0), 

for any (r, 0) e M+ x T and z e S T := {z e C : | Imz| < r}. 
(A5) v is differentiable with respect to r and K 9 a H> d r v(r,9 + a) can be 
analytically continued to d r v z {r,0) on C. Moreover, there exist a p > 
such that for any r > there is n T > such that \l^ r>p jd r v z (r, ff)\ < k t 
for any (r, fljel+xT and z£S T . 

Theorem 3. Assume that B satisfies (Al) and (A3) and V satisfies (A2),(A4) 
and (A5). Let A e if oc (R 2 ; R 2 ), p e (2,oo], with curlA = B. Then, for any 
eigenfunction ^> of H = Da + V with H^> = and E 6 M\ o~(Da ) the following 
holds: For every < a < 1, we have 

(6) ||e QBo/4|x|2 *|| 2 < oo. 

This theorem is proven in Section [5] 

Remark 5. The analyticity assumption (A4) on the angular variable ofV implies, 
by a Paley- Wiener argument, exponential decay of the Fourier modes of the potential 
in its angular momentum decomposition (see equations (|44l) . (1451) and (1611) below). 
The assumption (A5) is similar to (A4 ) but for the radial derivative of the potential. 

Remark 6. The first proof of Gaussian localization for magnetic Schrddinger op- 
erators using assumptions like (A4) (but not (A5)) was given in |10j . In addition, 
an example of a potential decaying at infinity for which the corresponding ground 
state decays slower than a Gaussian is also given in [10] . The proof in [10] is based 
on a generalized Feynman-Kac formula. An alternative proof using Agmon-type 
estimates with localizations in space and angular momentum was given in |19] . A 
variation of the method in |19j was used in |28j to treat the general n-dimensional 
case, again for magnetic Schrddinger operators. Our proof follows the ideas devel- 
oped in |19j . However, it turns out to be more involved since our operator is not 
bounded from below. To overcome this difficulty we square the Dirac operator ( or 
parts of it). This is the reason why (A5) is used in our setting. 

The article is organized as follows: In Section [3] we review some essentially well 
known facts about magnetic Dirac operators. Sections IH El and |6] are devoted to 
the proofs of theorems [T] [21 and [3] respectively. The article ends with an appendix 
containing some useful technical results. 

Acknowledgements. E.S. thanks Horia Cornean for stimulating discussions in the 
conference 'Spectral Days' in Santiago. Both authors have been partially supported 
by the DFG (SFB/TR12). 

3. Preliminaries 

3.1. Essential self-adjointness. Throughout this article we consider magnetic 
potentials A e if oc (R 2 ;M 2 ) and electric potentials V € L^ oc (R 2 ), p,q£ (2,oo]. In 



LOCALIZATION OF MASSLESS DIRAC FERMIONS 



5 



order to show essential self-adjointness of the operators H and Da defined in (p} 
and ([2]) it suffices to prove that 

H R <j> := a ■ (p - l { | x |< fl} A)0 + l {W < fl} ^, £ C °°(1R 2 ; C 2 ) , 

is essentially self-adjoint for every R > (see [5]). Using that for / G L P (R 2 ;C) 
and 2 < p < oo 

/(x)(p 2 + l)- 1 / 2 

is a compact operator (see [571 Theorem 4.1]) we get that l{| x |<_R}(l /r — <x • A) is a 
relative compact perturbation of Dq. This shows essential self-adjointness of Hr, 
since Do is essentially self-adjoint on C§°(R 2 ;C 2 ). 

3.2. Gauge invariance. Let A, A G Lf oc (R 2 ; M 2 ), 2 < p < oo, be two vector 
potentials with 

curl A = curl A 

in the sense of distributions. According to p~7j there is a gauge function $ G 
(M 2 ;R) such that 

A = A + V$. 

It follows, for any electric potential V G L^ oc (K 2 ; E), g € (2, oo], that 

(D A + y) = e i *(D A + U)e- i<l . 

In particular, e i *(D A + V^e" 1 * is essentially self-adjoint on C^°(R 2 ;C 2 ). 

This can be seen as follows: Note that V{e li {D k + V)e- it% ) = {/ G L 2 (M 2 ;C 2 ) : 

e -i$j g j){D A + V)}. Pick functions rj, rj G C£°(M 2 ; C 2 ) and a sequence ($ m ) meN 

in C°°(]R 2 ;IR) with $ m -> $ in W^(E. 2 ) (and hence in W^IR 2 )) as m -» oo. 

Then, 

((D A +yy ! e- i %)= lim ((D A + y) ?7 ',e- i *™7 7 ) 

= lim (e i *-7 ? ',(D A + - lim (e 1 <r • V$ m7? ) 

= (e 1 V, (D A + 7)r?) - (e 1 *rf, <x ■ V$r?> . 

Since ?/ is an arbitrary element of a core of D A + V, it follows that e -1 *7y G 
2?(D A + V) and (D A + V^e -1 *?? = er' lih {D k + V - a ■ V$)?? which implies that 

e l * (D A + V)e~ i *r? = (Da + V) V , t? g C%° (M 2 ; C 2 ) . 

Due to the essential self-adjointness of Da + U we deduce that e 1 *(D A + V r )e _1 * 
is also essentially self-adjoint on C^°(R 2 ;C 2 ) and that the two operators coincide. 

3.3. Supersymmetry. For A = (A 1 ,A 2 ) with Aj G Lf oc (K 2 ), p G (2, oo], j = 1, 2, 
we define the following two operators 

drf = [(p! - A x ) + i (pa - i4 a )]0, G C °°(]R 2 ; C) , 

d 2 ^> = [(p! - Ai) - i (pa - 4 a )]0, G C °°(]R 2 ; C) . 

Clearly, we have that 

Da rc -(H*i(?)= ( I d Q J ■ 
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Since Da \c°°(r 2 ;C 2 ) is essentially self-adjoint it follows that d\ and di are closable 
[2"9l Section 5.2.2]. In addition, setting d := d\ one finds that d* = d 2 and 

(7) DA= (d 0*) ° n v {D A )=V{d)®V{d*). 

It is known that dd* and d*d are self-adjoint with domains T>(dd*) ={(/)€ T)(d*) : 
d*4> e V(d)} and V(d*d) = {4> € 23(d) : d<f> e X>(d* )}. Moreover, there is a unitary 
map S from Ker(dd*) ± to Ker(d*d)- L , such that 

(8) dd* |'Kcr(dd*)- L = S*d*d tKcrCd'd)- 1 - ^ ■ 

Let us note that we can block-diagonalize Da using the Foldy-Wouthuysen trans- 
formation. Setting 

1/V2 on Ker^A)- 1 f l/\/2 on Ker(Z?A)- L 

1 on Ker(£> A ) ' °~ ~~ \ on Ker(L> A ) ' 

we define the Foldy-Wouthuysen transformation as 

U = a + + <7 3 sgn(D A )a_ , 
where sgn(Z?A) = Da/\Da\ on Ker(Z?A) ± and equals zero on Ker(£>A) and 

1 



"'■'■ = ^ -1 

The unitarity of the above transformation can be easily verified observing that 
(p E Ker(DA) o~3<P 6 Ker(£>A) and that a^sgn(DA) = — sgn(Z?A) cr 3- The latter 
relation holds since o^Da = -Ba^ and CT3|Da| = |-Da|<73. A direct computation 
yields 

fdFd 
-\fdd* 



(9) UDaU* = 
Equation © and ([§]) imply the following statement. 

Proposition 1. Let A € Lf oc (R 2 ;R 2 ) for some p £ (2,oo]. Then, the spectrum of 
Da is symmetric with respect to zero and 

o-j {Da) n (0, 00) = (Tj) ( VcFd) \ {0} , fl e {pp, d, ess} . 
4. The spectrum of Z? a 

The aim of this section is to show Theorem [TJ An important ingredient is the 
study of the essential spectrum of D A . In order to do that we modify an argument 
from [TJ] obtaining Lemma [JJ below. We combine this with a result from [53] on 
the infiniteness of zero modes for Pauli-operators (see Lemma [5] below) . The proof 
of the theorem is then a consequence of lemmas and [31 

In the following discussion we assume that B = Bq + b with Bq > and b E 
L 1 1 oc (M 2 ; M) such that {bl 1 / 2 is relative a/p 2 + 1- compact. Let A € Lf oc (R 2 ; R 2 ), p £ 
(2, 00], with curl A = B. We start by observing that, for cj> e Cg° (R 2 ; C), 

2 

(10) j t 

2 



= £ \\( Pj -A. 
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holds. This implies the commutator relation 

(11) {(j>,[d,d*](t)):={<f(j>,d*<j))-{d^d(t))=2{^B(t)), e C^{R 2 ; C) . 

The idea in |14] is to use this commutator to study the essential spectrum of dd* 
and d*d. In order to extend this identity we define these operators as quadratic 
forms and show that Q(b) D Q(d*d) = Q(dd*) and |fe| 1/2 (d*d + l)" 1 / 2 is a compact 
operator. Here Q(-) is used to denote the form domain. 
Let us define 

qi {<t>,4>) =q 1 [<f>] := ||#f , g 2 (<^) = ~ ||d>|| 2 , 

with form domains Q{qi) = 13(d) and Q(q2) = T>(d*). Since d and d* are closed (see 
Subsection l3.3[) we have that qi and q2 are closed and positive. Thus, associated to 
qj, j = 1, 2, there is a unique self-adjoint operator Tj characterized as follows: 

' ' 2 ' 2?(T 3 -) = e Qfe) | 377 e L 2 (R 2 ; C), W> e C, fcty, = (V, v)} , 

where C is any form core of qj. It is easy to check using (jT2"j) that in fact T\ = d*d 
and T2 = dd* . Note that since the restrictions of d and d* to (R 2 , C) are closable 
C5°(M 2 , C) is a form core for (71 and 92. We define yet another quadratic form. For 

4> € Cq° (R 2 , C) we set 

2 

It is known [26] that 93 is closable and we denote its closure by g 3 . Its associated 
self-adjoint operator H$ =: (p — A) 2 is the usual magnetic Schrddinger operator. 
Recall that |6| 1//2 is relative \/p 2 + 1-compact. Using the diamagnetic inequality 
for |p — A| (see e.g. [11]) and arguing as in [2j Theorem 2.6] we conclude that 
Q(q 3 ) = V(Hy 2 ) C £>(|&| 1/2 ) and that I&I 1 / 2 is relative id^-compact. Thus, the 
quadratic form 

/%] :=i?o||^|| 2 + (sgn(fe)|6| 1 /V,|6| 1 / 2 ^) 
is in absolute value bounded with respect to 93 with bound 0. In particular, 

<if[v} '■= 93 M ± /%] ) peSfe), 

is closed. Observing that by (|10[) we have that 92 \c° a= 9^ \c°° and 91 tc°° = 
93~ \cg° and using that Cg°(R 2 ;C) is a form core for (71,92,93 and g 3 we conclude 
that Q(qi) = 2(92) = 2(93) = Q and 91 = 9^ and 92 = 9^. Moreover, 

dd* = (p — A) 2 + B 

(13) 

d*d= (p- A) 2 -£ 

in the sense of quadratic forms on Q and hence the commutator formula (jlip 
extends to Q. 

Lemma 1. Let B = B + b with B > and I6I 1 / 2 G L 2 oc (R 2 ;R) 6e retowe 
V^P 2 + 1-compact. Let A <G Lf oc (R 2 ; R 2 ),p € (2, 00], udi/i curl A = B. Then, either 
one of the following statements holds 

i) a ess (d*d) = 

ii) cr ess (d*d) = {2B n : n e N } and a ess (dd*) = {2B n : n G N} . 



s 



MARTIN KONENBERG AND EDGARDO STOCKMEYER 



In addition, if V satisfies (A2) then V is relative Da- compact and in particular 
CTcss(-Da) = cr C ss (-Da + V). 

Remark 7. Note that our assumption on B are satisfied if B fulfills (Al). Indeed, 
in this case |6| 1/,2 (p 2 + l) -1 / 2 is compact by Lemma \T2A in Appendix ["Si 

Moreover, note that if B G Lj 9 oc (R 2 ; R) for some q > 1 we can always find 
A G -L^ oc (R 2 ; R 2 ) for some p G (2, oo]. In order to see this define h to be a solution 
of 

(14) Ah = B. 

A local solution to this equation is given by the Newton potential hjy of B. We 
know that hfj G W 2,q (fl) by the C alder on- Zygmund inequality, where O C R 2 is a 
bounded domain (see e.g. p21 Sec. 9. 4]). This property extends to any solution h of 
(|14p since h — h,N is harmonic on fl. Therefore, h G W,'£(M, 2 ). Now one can define 
A := (— c^/i, Clearly, Aj G W„'5(K 2 ). By standard Sobolev inequalities one 

obtains that A j G L* oc (R 2 ;R) for some 2 < t < oo if q G (1, 2] and A; G L£ C (R 2 ;M) 
i/g>2. 

Proof. First note that for any A > the operator (d*d+2B +\y 1/2 maps L 2 (R 2 ; C) 

1 /2 

onto T>(vd*d) which equals Q and T>(H S ' ). Thus, by the closed graph theorem, 
the operator (H s + l) 1 / 2 (d*d + 2B + A) -1 / 2 is bounded. In particular, 

\b\ l ' 2 {d*d + 2B + X)- 1 ' 2 = |6| 1/2 (i/ s + iy 1/2 {H s + l) 1 / 2 {d*d + 2B + A)" 1 / 2 

is compact. Hence, the operator 

T(A) := (d*d + 2B + X)- 1 / 2 sgn(b)\b\ 1 / 2 \b\ 1 / 2 (d*d + 2B + A)^ 1/2 

is also compact. It is easy to see that A > can be chosen so large that ||T(A)j| < 1. 
For such A's we have, according to the resolvent formula for operators defined as 
quadratic forms (see [25]), that 

(d*d + 2B a + 2b + X)- 1 = (d*d + 2B Q + A)- 1/2 (l + T(A))" 1 (d*d + 2B Q + A)~ 1/2 . 

Note that the inverse of 1 + T(X) is well defined as a geometric expansion. Since 
(1 + T(A)) -1 — 1 is compact, we conclude that the resolvent difference between 
d*d + 2Bq + 2b + X and d*d + 2Bq + A is also compact. Therefore, by Weyl's 
theorem, the two operators have the same essential spectrum. Using this and (fT5|) 
we deduce that 

(15) a css (dd*) = a css (d*d + 2B + 2b) = a css (d*d + 2B Q ) . 
The latter equality and Equation (|5J) imply (here we follow |14j ) 

S := a css (d*d) , S C [0, oo) , 

(16) S\{0} = S + 2B Q . 

Assume now that S ^ 0, then it is easy to see from (fit?]) that G S and hence 
2B n G S, n G N . Note also that no other points can belong to S. Hence, using 
(fl"5)) we get that a css (dd*) = 2B n, n G N. 

Now, assume that V fulfills (A2). Then, V is relative \/p 2 + 1- compact (see 
Lemma [12] in Appendix |B]) . It follows by the diamagnetic inequality that V is 

1 /2 

relative H s ' - compact and consequently (arguing as before for b) V(V) D Q and 
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the operators V(dd* + A 2 ) 1 / 2 and V(d*d + A 2 ) x l 2 are compact for any A ^ 0. 
From these considerations follow that V(D A — i A) -1 is compact, since the identity 

(D A i A)- 1 = (D% + X 2 )- 1/2 [(D 2 A + X 2 )- 1/2 (Da + i A)] 

/ (d*d + X 2 )- 1 / 2 \ 2 ,^-l/2 (n I : \\1 

-[ o {dd* + \ 2 r^ ) -K^a + a ) (-Ca + A)] 

holds and the operator in [. . . ] is bounded. Therefore, 

Vess{D A + V) = a ess (D A ). 

□ 

We note that if b satisfies (Al) then Ker(d*d) is infinitely degenerated. Indeed, 
this follows from the fact that 

(17) / [B) + d 2 x = oo, / [B]_d 2 x <oo, 
Jr 2 Jr 2 

(where [/] + and [/]_ are the positive and negative parts of /) which shows that 
B = Bq + b fulfills the conditions of [531 Corollary 3.4]. In particular, we know that 

Ker(d*d) = {uie~ h | we - ' 1 € L 2 (R 2 ; C) , cj is analytic in x\ + ix 2 } , 

where h is a solution of the equation Ah = B |24| . Therefore, we get: 

Lemma 2. Assume that B satisfies (Al) and let A <E Lf oc (]R 2 ; R 2 ), p <E (2,oo], 
with curl A = B . Then, 

(18) a css {d*d) = {2S n|ne N } , a css (dd*) — {2B n | n G N} . 
In particular, 

cr css (D A ) = a css (D Aa ) = {l n | n € Z} . 
Moreover, is an isolated point of a{D A ) and a(d*d). 

Proof. Due to our previous discussion we see that G a ess {d*d). This combined 
with Lemma Q] imply (fT8]) . That is an isolated point of o~(d*d) follows by noting 
that, since a ess [dd*), is neither an accumulation point of a(dd*) nor of a(d*d). 
The statements on cr(D A ) are now a consequence of Proposition [TJ □ 

Lemma 3. Assume that B satisfies (Al) and let A e i^ oc (]R 2 ; R 2 ), p e (2, oo], 
with curl A = B. Then, we have: 

(a) If b < and strictly negative on some open set, then 

dim(Ran(l (0V ^ ) (D A ))) = dim(Ran(l ( _ V 2 B?i0) (Z> A ))) = oo . 

(b) lfb>0 then 

dim(Ran(l (0v ^ ) (I? A ))) = dim(Ran(l ( _ v ^ 0) (C A ))) =0. 

Proof. We may choose A := {—d^h, d\h) where h is a solution of Ah = B. Due to 
Remark[7]we know that A e Lf oc (R 2 ,K 2 ) for some p > 2. Part (a): Let il be an 
open set with b \ £1 < 0. Recall that there are infinitely many functions uj, analytic 
in xi + 1x2, with ip := use~ h 6 Kcr(d*d). For such ip we have, using (fT3|) . 

(19) dd» = 2<^>, Bil>) < 2S ||V'ir + 2 / 6(x)|V>(x)| 2 dx < 2S || V|| 2 , 

Jo. 
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where in the last inequality we use the fact that ip can not vanish on Q. Let (V'n)neN 
be an orthonormal system such that ip n '■= er h uj n G Ker d* d with Lu n analytic in 
xi + ix2- For N G N define the self-adjoint matrix Mn '■= ((f/Vt, dd*ip m ))i< ntm <N ■ 
It follows from (fl9|) that Mat < 2Bq. The Rayleigh-Ritz principle implies 

< Hn{d<f) < »n(M N ) < 2B , n = 1, . . . , N , 

where we write 

H n {T) := sup inf (ip,Tip) 

<fll,...,(f>n — l *S s P an {*l''-'>*n-l}' 
||*||=l,</>6QCr) 

for some self-adjoint operator T. Since N is arbitrary the mini-max principle implies 
that dim(Ran(lr y /2B^){dd*))) — oo. It follows that dim(Ran(lj ^2B^)(dd*))) — 
oo, for ^ a css (dd*) by Lemma O The claim is now a consequence of Proposition 
□ and ©. 

Part (b): In this case we have that dd* > 2B , since dd* — d*d = 2B > 2B . Thus, 
the claim follows now from Proposition Q] and ©. □ 

5. Super-exponential localization 

The proof of Theorcm[2]follows the ideas developed in [6]. An essential ingredient 
is that, by means of suitable local gauge transformations on certain regions outside 
a big ball of radius n centered at the origin, one can replace the operator D A by a 
Dirac operator Z?a„ with A n = Ao + a™ , where a™ is a magnetic vector potential 
of a magnetic field b n satisfying linin^oo ||6 n ||oo = 0. The advantage of this is 
that we can obtain explicit L p estimates (see Lemma [4] below) for the resolvents of 
Da„ i conjugated with exponential weights. These estimates can be derived using 
a certain resolvent expansion (see (|38p ) in combination with an explicit expression 
for the Green kernel of £>a that can be found in Appendix [X] below. 

Before stating these LP estimates let us fix some notation. For p, q € [l,oo] we 
denote by B(p,q) the space of bounded operators from L P (R 2 ;C 2 ) to L q (R 2 ;C 2 ) 
and write, for T € B(p, q), 

(20) ||T|| p , g := \\T\\ BM . 

Let 7 > and u€E 2 with |u| = 1. We define the exponential weight function as 

F(x):= 7 u-x, xel 2 . 

Let b n be a magnetic field with lmin-^oo ||6 n ||oo = and a™ be the associated vector 
potential in the transversal gauge, i.e., 

(21) a n (x) := / 6„(sx)Axsds, 

where we write a A v := a(— V2, vi) for a € R and v S R 2 . The proof of the Lemma 
below can be found at the end of this section. 

Lemma 4. Let V n G L°°(R 2 ; R), n G N, be a family of electric potentials satisfying 

lim HKlIco =0. 

n— >oo 

For any n G N define the family of self-adjoint operators Z?a„ + V n , where A n := 
Ao + a" and a" is given in (|21[) . Let z G R \ ct(Da„) and q,r € [l,oo] be such 
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that 1 + p — i = | /or some p G [1, 2). Then, there exists N > suc/i i/iaf, /or aZZ 
n> N , z ^ u(Da„ + K) 

(22) e F ( J DA„+K-^" 1 e- F eB( g ,r). 

In what follows we apply the above result to show Theorem [5] 
Proof of Theorem^ For n G N and u G R 2 with |u| = 1 set 

(23) a n = {x G R 2 : u • x > n} . 

For j G {1, 2, 3} define Xj € C°°(R 2 ; [0, 1]) with Xj = on R 2 \ fl jn and Xj = 1 on 
Qrj + X) n . We choose n so large that 

IHU°°(n n ) < 00 • 

Since 6 € L^R 2 ), q > 1, we find a vector potential a G Lf oc (R 2 ;R 2 ), p > 2, with 
curia = 6 (see Remark [7]). Define, for x G R 2 , 

a n (x) = / b n (sx) A xl sds , 
Jo 

where b n :— ln n b G L°°(R 2 ). Observe that 

curl a = curl a" on f2„ , 

that f2 n is simply connected, and that a", a G L^ oc (R 2 ; R 2 ) for some p > 2. There- 
fore, there exists a gauge function $„ G W lo ',?(fi„) such that (see [T71 Lemma 1.1]) 

(24) V$„ = a a" on ft„ . 

By multiplying $„ with a C°°- cutoff function we may define a $„ G Wf'^ (R 2 ) that 
coincides with $ n on f^n- In particular, we find that 

(25) V$„ = a a" on n 2n • 

Define now V n :— XiV and observe that ||14,||oo — > as n — > oo. Then we get, for 
any 77 G C^°(R 2 ; C 2 ) and j G {2, 3}, using (f2"5"j) and the identity Xj = XiXj> 

Xi(-°A + V)T} = X 3 {Da„ er • a + K)?7 

= Xj(-C»Ao - ■ V$„ - <t ■ a" + V^)»j 
= X 3 e [ ^(D Ao ~ a ■ a n + V n )e-' 1 ^r 1 . 

Set A n := Ao + a™ and let ^ be an eigenfunction of D A with eigenvalue E ^ 
dpAo)- By the previous computation we obtain, for any r\ G Co°(R 2 ;C 2 ) and 

je{2,3}, " 

( e l (£> An + K - £)e- ! *" ry | Xj * > = ( (D A + V - E) V \ Xj * ) 

= (i<r- Vx^l*). 

This equality extends to any r\ in the domain of e'*"(-D Ajl + Ki — £ , )e -1 *", since 
Cg°(R 2 ;C 2 ) is a core for e i<E> "(L> A „ + K - E)e~ i ^ (see Subsection EH} . Clearly, 
A„ and satisfy the assumptions of Lemma £f] Thus, for n sufficiently large, 
E ^ cr(D Atl ) and we may replace 77 by e 1 * n (-D A „ + V n — E)~ 1 e~ 1 *™?7 obtaining that 

(26) x J * = -ie i *"(^ Al +K-^)" 1 e M *"(T-Vx 3 )*, jG {2,3}. 
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Observe that using (|26|) for j = 2 in combination with Lemma |4] (with q = 2, r = 3, 
and F = 0) we obtain that 

(27) X2* € L 3 (R 2 ;C 2 ). 
We use again (|26p, for large n, to get in addition that 

(28) e f X 3* = -ie 1 *" (e F (£> An + K - i?)" 1 ^) ( e - 1 *»e F (<rVx3)X2*) • 



Since supp(Vxs) C ^3 n \ £li n we find thanks to ([2?) that e 1 *"e i? (crVx3)x2* G 
L 2 (R 2 ;C 2 ) n L 3 (R 2 ;C 2 ). Thus, we may apply Lemma H with q = 3,r = oo and 
q = 2, r = 2 to obtain the decay in the L°° and L 2 norms respectively for n > uq 
sufficiently large. We obtain the desired bound (|5|) from (f2"5| by varying F over 
finitely many vectors u. □ 

Proof of Lemma [^J Recall that the magnetic vector potential is given by A„ = 
Ao + a™ where a" is defined in (|21|) . 

A simple calculation shows that the vector potential 



(29) 



a x'( x )= / b n (x.' + s(x — x')) A (x — x') sds , x' £ 



is also a vector potential of the magnetic field b n . A crucial property of a™, is that 

(30) |b£,(x)| < II^IU-lx-x'l, x,x'eR 2 . 
Since curl a£, = curl a 11 there exists a function (p n : R 2 — > R with 

(31) V x ^ n (x,x')=a"(x)-a^,(x). 
We may further require that 

(32) ¥ >„(x,x) = 0. 

The proof of Lemma |4] is based upon L p estimates for the resolvent expansion (|38|) 
below. We start by defining the relevant objects and list their LP properties. For 
z £ R\cr(£>A ) let Go(x, x', z) be a representation of the Green kernel of (Da — z)~ x 
as 2 x 2-matrix. Due to (|81l) from Appendix [X] and the triangular inequality we 
obtain that 

\\e F ^G (x, x'; z)e- p ^ || cw < e^^+^-'L^ ~ *i z) := 5 f(x - x') . 



We observe that thanks to (|82|) we have that 

(33) g F G L'(l 2 ), i€[l,2) and \x\g F £ L^R 2 ), t£[l,oo]. 

We introduce, for n £ N, the integral operators S n (z),T n (z) : L 2 (R 2 ;C 2 ) 
L 2 (R 2 ;C 2 ) with 

(34) (S n (z)/)(x):= / e i ^^ x ')G (x,x';z)/(x')dx', 

JR 2 

(35) (T„(z)/)(x):= / < T.a^(x) e i ^^ x ')G (x,x';z)/(x')dx', 



where y>„ is determined by (|5Tj) and (j3"2")l . Notice that in view of (f3"3")l . (1501) . 
and Young's inequality (see [HI Section 4.2]) both operators are well defined and 
bounded. In fact, since 

||(e F T n (z)e -i? )(x,x')||c3®C= < ||6„|| oc |x - x'|5 F (x - x') , 
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we find by (|33|) and Young's inequality that, for q £ [l,oo], 
(36) Km ||e F T n (z)e- F ||,, 9 = 0. 



Furthermore, a similar argument implies that, for q,r £ [l,oo] and t £ [1,2) with 
1=1+1-1 

(37) sup ||e- F S'„(z)e _F || 9ir < oo . 

Our next task is to show the following resolvent formula in L 2 (R 2 ; C 2 ), for n E N 
so large that ||T„(z)|| 2 ,2 < 1 (see ([36])). 

OO 

(38) (D An - z)- 1 = S n (z)J2T n (z) k . 

k=0 

Pick functions / G L 2 (R 2 ; C 2 ) and g G C^(M 2 ; C 2 ) and, we find 
((D An -z)g,S n (z)f) 

= [ ([(D An - z)g](x), [ e i ^^ x ')Go(x,x'; Z )/(x')dx') c2 dx 
Jm 2 Jk 2 

( e -i<Mx,x') [(_Da„ - z) 5 ](x), G (x,x'; z)/(x')) c2 dxdx' , 

where Young's inequality together with Lemma QT] from Appendix 1X1 enabled us to 
use Fubini's theorem in the last equality. Observe that due to (f5T|) 

e - ly „(x,x') rp^ _ z)5](x) = p^ _ ffia j_ z)e -i(P-(-.x') ff](x) . 

Hence, using again Fubini's theorem, 
((D An -z)g,S n (z)f) = 

( f Go(x',x;z)[p Ao -z)e- i ^(- x '). 9 ](x)dx,/(x')) c2 dx' 

(<r • a" (x)e- i ^^ x ').g(x), G (x, x'; z)/(x')) c2 dx' dx 

= (9,f)-(g,T n (z)f). 

For the first integral after the first equality above we used (I52"j) together with the 
fact that Go is the Green kernel of D Ao and thus, for any g £ G 3O (M 2 ; C 2 ), 

/ Go(x',x;z)[(£) Ao - ^)5]( x ) <^ x = 5( x ') a.e. . 

Jm. 2 

Now, since -Da„ is essentially self-adjoint on G X> (1R 2 ; C 2 ) we can extend the identity 
(1551) for all 5 € V(D An ). From this extension follows that maps L 2 (R 2 ;C 2 ) 

in T>(D An ) and 

p A „ - = / - T„(z)/ , / e L 2 (R 2 ; C 2 ) . 

This yields, for n large enough, the operator identity 

S n (z) = (D An -z)- 1 (l-T n (z)), 

from which follows 



(39) 
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We observe that (|22|) . for V n = 0, is a consequence of (|36|) . (|37|) . and 

oo , 

e F (D An -z)-'e- F < e F S n (z)e- F ■ V (\\e F T n (z)e- F |l ) 



q.r 

H ' k=0 



Note that the last sum above converges for n large enough due to (|36D . 
In order to show (|22|) for V„ =^ we note that by Holder's inequality 

(40) \\V n e F {D An - z)- 1 e- F \\ q , q <\\V n \U \\e F (D An - z)- x e- F \\ q>q -► , 

as ri — >■ oo. Therefore, the following computation is meaningful for n large enough 

.fm. , t/ „\-i„--Fi 



K (Da, + V n - z)- L e 



I 9)^ 



= \\e F (D An - z)-\l + V n (D An - z)- 1 )- 1 er F \\_ 

OO 

< \\e F (D An -z)- 1 e- F \\ q , r \\{V n e F (D An -z)- 1 e- F } m \\ q , q . 
This finishes the proof. □ 



6. Gaussian-localization 

In this section we show Theorem [3] on Gaussian localization of eigenfunctions 
with energies in the discrete spectrum of 

(41) H = D A + V, 

under the assumptions (Al)-(A5) stated in the introduction. We choose the mag- 
netic vector potential to be given by 

(42) A(x) := r^Air) ( ^f 2 ^ , A(r) = r" 1 J B{s)s ds . 

If B E L|' oc (K 2 ,M) it is easy to see, using Holder's inequality, that if q € (1,2] 
then A e Lf oc (M 2 ; M 2 ), for some p € (2,oo), and that A G L^ C (M 2 ;R 2 ) whenever 
q € (2,oo] . 

The proof of Theorem [31 given in Subsection 16.31 follows the ideas of QI5] con- 
sisting in Agmon-type estimates with localizations in space and in the angular 
momentum variable. Of course, we have to adapt the method of |19j since our 
Hamiltonian is not bounded form below. 

In Subsection 16.11 we transform the operator H to polar coordinates and we 
decompose it in the angular momentum variable rrij. The analyticity condition 
(A4) on V permits us to obtain exponential decay in \rrij\ of eigenfunctions of H 
with eigenvalues E £ <Jd{H) (see Lemma[5]in Subsection l6.2l) . In order to obtain the 
Agmon estimates, in Subsection l6.3[ we square the transformed free Dirac operator 
Kq (see (j43|) for its definition) . The Gaussian decay is essentially due to a positive 
term in (Kq 2 ^) 2 that goes like r 2 . This term is in competition with a term that 
behaves like rrij when rrij > 0. The Gaussian weights in the Agmon estimates are 
localized in the region where rrij < r 2 . The complementary region, on the other 
hand, is controlled by the exponential decay in \rrij\. 
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6.1. Unitary transform. In the following we derive an equivalent representation 
of H . We denote by U the unitary map that represents H in polar coordinates (see 
e.g. [25J Section 7.3.3]) 

UHU* = H {1) = + v(r, 9), 

:=Se{-id r + ir- 1 a 3 Jz-Y ( T 5 A{r)}, 



acting on U {1) := L 2 (R+) ® L 2 (T; C 2 ) 2 , where 
J 3 := -ide + l/2o-3, S e := 



e" i( 
s ie 



Next we identify L 2 {T; C 2 ) with l 2 (Z) 2 by means of the transformation T : L 2 (T; C 2 ) 
£ 2 (Z) 2 given by 



_ r 

2?r 



: = -L / Mge- im ^f(6)d9, 



for / G L 2 (T; d6>) 2 , where mj = (2j + l)/2, j e Z, and 

M - ( ei9/2 

Under these transformations we find the decomposition 

L 2 (M 2 ; C 2 ) = H (2 > := L 2 (R+; dr) 2 

and the corresponding operator 

H = H (2 > = i^ 2) + tU , 

which is essentially self-adjoint on := J"J7C^°(R 2 ; C 2 ). For /i G X>( 2 \ i^ 2) = 
TUD A U*T* acts as 

(43) (K^h)(r,j) = ( - i<j 2 r + a^m.r- 1 + A(r)))h(r,j) , 

where we used that TSeJ 7 * = &2, TSga^J-* = i<7i, and that TJ3T* is the multi- 
plication operator by rrij. The electric potential W = FvF* acts as 

(44) (Wh)(r,l) :=Y,i>(. r > l -3)Kr,j), 

jew. 

where 



(45) v(r,n) = — e - ine v{r,9)d9 , neZ. 

27r Jo 

Two other quantities play an important role in our analysis, namely W\ '■= J-d r vJ-* 
given by 

(46) (Wih)(r,l) :=J2 d rt(r,l-j)h(r,j), 

jez 

and W2 ■= J-devF* that acts as 

(47) (W 2 h)(r,l) :=J2i(i-l)HrJ-3)h(r,j). 

jei 
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6.2. Rotation-analyticity. For / £ and a£lwe set 

(48) (U a f)(r,9) := (e ij3 °/)(r, 0) = j°* a ' 2 f(r,6 + a) . 
We call a vector / £ T^ 1 ) rotation-analytic, if and only if the series 

' n! 

has an infinite radius of convergence. We start by presenting a lemma that gives 
us some a priori decay of some eigenfunctions of H ^ in the angular momentum 
variable. 

Lemma 5. Assume that (Al)-(A^) hold. Let £ be an eigenfunction of 
to the eigenvalue E £ ad(H^). Then, for every 7 > 0, we have 

P 00 

(49) V / e 2 ^ m ^\*(r,j)\ 2 dr < 00. 

Proof. The proof is analog to the one given in [TH Section 3] . We sketch it here for 
the reader's convenience. Due to Lemma fl4l (in Appendix |B|) {H^(z)} ze c defined 
on V(K { 1} ) through 

is an analytic family of type (A) (see [13]). Note that when a £ K the identity 
H^\a) = U a H^U* holds. Moreover, by Lemma fT3l (in AppendixEJ we have that 

is a compact operator in for any z £ C Therefore, a ess (H W (z)) — cj css (.K"q 1 ' ) ) 
by Weyl's theorem. Arguing with analytic perturbation theory and using the fact 
that the spectrum of ifW(a) and is the same for a real (see e.g. the proof 
of Theorem XIII. 36 in 23] for a similar argument) we find that E £ a^H^) of 
multiplicity N £ N is also an eigenvalue of H W (z) of the same multiplicity. 

Let P z be the iV-dimensional P-eigenprojection of H^(z). Since rotation- 
analytic vectors are dense in (see e.g. [20]) we find some rotation-analytic 
vectors /1, . . . , /jv such that RanP = Span{P /i, . . . , -Po/jv} ■ Observing that, for 
a £ R and j £ {1, . . . , N}, 

U a Pofj = P a U a fj , 

we find an analytic continuation of a h> U a Pofj £ to the complex plane. In 
particular, e lj3Z Pofj belongs to US 1 ' for any z £ C. Let € RanPo be such that 
J^i = ^P. By the discussion above we get that 

This ends the proof since (J r e j37 \E , i)(r, j) = e" lj ^^(r, j) and 
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6.3. Agmon-type Estimate. In this section we deduce the Agmon estimates 
needed in the proof of Theorem [3] They were obtained in [19] for the case of 
a magnetic Schrodinger operator. These estimates uses heavily the fact that the 
Schrodinger operator is bounded from below. As we commented before we will 
obtain these estimates for the square of the Dirac operator Kg . 

Fix a number B > Bo and note that, due to (A2), there exists i?o > so large 
that the estimate (l58l) is fulfilled and moreover 



(50) \\l {r>Ro} B\\ < B , r > Rq . 

We set, for < (fe < Qx < 1, 



(51) r(j) ;={ V(A^'' m J-° 

[0, nij < 0, 

(52) n quq2 :={(r,i)GK+xZ|r>r(i)}. 
Moreover, we define 

(q 2 B /4(r 2 - r(j) 2 ), m 3 > 0, r > r(j) 

(53) P(r,j) := lq 2 B r 2 /4: mj < 

[0 rrij > 0, r < r(j) . 

Eventually we will choose q 2 to be sufficiently close to 1. A direct calculation shows 
that 



(54) \p(r,h)- P (r,j 2 )\ < - 2 ^ 2 ,„ \ji - j 2 \ 



q 2 B 
(q 2 i-q 2 2 )Br, 

Let p e := p(l + £p)~ l • It is easy to see that 



(55) \pe(r,jx) -pe(r,h)\ < . 2 2 2 , \h~h\- 

[q x - q 2 )B Q 

Finally, for R > i?o, we fix a smooth function in r with bounded derivatives in 
R+ x Z satisfying 



(56) f R (r,j) 



1 r > 2R and (r, j) G 0. quq% 
r<i? or (r,j) gn quM2 



where [i € (0, 1) is a fixed number that will be chosen sufficiently close to 1. Note 

that Qqi,q2 ^Qi,/^92' 

Let \& be the eigenfunction from Theorem [3] and $ := J-U^f be a normalized 
eigenfunction of H^> with corresponding energy E e a d (H^). We set, for R> R 
and <5 G (0, ^i), 

(57) 3 :=e^/ fl *. 
Observe that 6 can be chosen arbitrarily close to 1. 

Lemma 6. We find constants R\ > Rq and c > such that, for any 6 G (—1, 1), p > 
Ri, and j € {0, 1, 2}, 

(58) sup\\6 p e s pc W je ~ s P '9 P \\ <c, 

e>0 
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where Wo ■= W and 9 P := l{ r>p }. In particular, the commutator 

(59) [K { 2 \W Q ] = -i<j 2 W 1 + l -^-W 2 
satisfies the estimate 

(60) supHe^Mi^WoJe-^i <2c, R > Ri . 

£>0 

Proof. We show (|58|) only for j = 2, for the other cases follow analogously. For any 
m € Z and r > we define the analytic function C9z4 h m (r, z) := e~ lmz v z (r, 0). 
Using (|45[) and the decay and analyticity assumptions on i> (A5) we find for any 
r > p (sufficiently large), m G Z, and 7 € E that there is a constant C > such 
that 



1 /-27T ^ '' 27r 



(61) 



\v(r,m)\=—\ h m (r,6)de\ = —\ I h m (r, 9 - ij)d0\ 



Here we also used Cauchy's integral theorem and the fact that h m (r,z) is 2ir- 
periodic with respect to Re(z). In particular, replacing 7 by —7 in the above 
estimate we see that, for 7 > and m <E Z, the bound 9 p \v(r,m)\ < Ce~l m l 7 holds. 
Therefore, using ([55]) . (|47|) . and Young's inequality for £ 2 (Z 2 ;C 2 ) in combination 
with the Cauchy-Schwarz inequality for L 2 ((0, 00); C), we get, for 7 sufficiently 
large and every / G H {2 \ that |(/, 9 p e 5p < W 2 e~ 5 p '9 p f)\ is bounded by 



rOO 

/ E \f( r > 01 E Mfr * - J')H* - i|e' T|, - J ' l |/(r > i)|dr < C||/|| 2 



for some constant C > 0, where r := , a 92 f, D . 

Equation ([6"U| follows from ([55| and (j5"9")) . Equation (j59")l is a consequence of 

(62) [X< 2) ,H/]=J-[^ 1) ,«]J-* 

(63) =F(-i Sgd r v + ^-d v)r* , 

r 

and the fact that J-SgJ-* — o~ 2 and J-S0O-3J-* —ia\. □ 

Before continuing let us state a simple technical result. 
Lemma 7. For any 7 G E we have that e 7 ^/ fi * G V(K^ 2) ). 

Proof. Let A > and 77 G .Ft/C^IR 2 ; C 2 ). First observe that a simple computation 
shows that 

{d r e-< p 'f R )e- Xr 

extends to a bounded operator on %( 2 \ In addition, e Xr ^ G 'H^ 2 -' by Theorem [2] 
Therefore, we get by explicit calculation on FUC^ 2 (M 2 ; C 2 ), that 
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Since r\ can be chosen arbitrarily from the domain of essential self-adjointness of 
Kq ' we get the desired result. □ 

An important role in our analysis is played by the quantity 

(64) Q:=Rc(l4 2) e 5 ^g\K^e- s ^g), 

which is well defined due to Lemma [7] Before we show Theorem [3] we state two 
preparatory lemmata whose proofs are given in the next subsection. 

Lemma 8. There are R, e -independent constants C\,C% > such that, for R > R\ 
sufficiently large, 

(65) Q>(C 1 R 2 -C 2 )\\g\\ 2 . 

Lemma 9. There is an R, e-independent constant C3 and an e-independent con- 
stant C(R) such that, for R > R\ sufficiently large, 

(66) Q<C 3 \\g\\ 2 + C(R)\\g\\. 

Proof of Theorem^ Fix S, 51,52 S (0, 1)- Combining Lemma H] and [5] we find, for 
R > i?i sufficiently large, 

(67) ||.g|| < (dR 2 - C 2 - C^CiR) 

Since the right hand side of (|67p is independent of e we obtain, by the monotone 
convergence theorem, 

(68) \\e s "n 2 = lim || e ^*|| 2 < (sup ||.g|| + ||e 4 '(l - f R )\\) 2 < 00 . 
For M > 1 define 

(69) ^Wm = {(r, 3) G K + x Z I r 2 > Mr(j) 2 } . 

We have that £lq lt q 3 ,M C Q qi ,q 3 - Thus, for any (r, j) € Cl qi>g2! M, we get 

(70) p(rJ)= q -^(r 2 -r( j f)> q -f L i^^Y- 
Therefore, setting a := ^52(1 — M^ 1 ), we obtain 

(71) \\e aBo/4r2 l^ ill < co. 

If (r, j) £ £l qi , q2 ,M then 

(72) mj >^-® B » r2 =:[3r\ 

J ~ 4MB 

Thus, thanks to Lemma [5] we deduce, for any 7 > 0, that 



(73) \\e^ r l fi0 ¥||<oo. 



Choosing7 = a/f3-Bo/4 and combining (|73l) with (J7TJ) we conclude that \\ e aB °/ ir *|| < 
00. The latter holds for a > arbitrarily close to 1, since 5 and 52 can be chosen 
arbitrarily close to 1 and M > 1 can be as large as we want. This proves the 
theorem. □ 
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6.4. Proof of lemmas [8] and Q51 Before we give the proof of lemmas |8] and [9] we 

need a preparatory result. 

Lemma 10. For R > R\ sufficiently large we have that 

(74) \\K^g\\ 2 W2) > ^q 2 2 B 2 \\rg\\ 2 /4 - \\r^ g\\ 2 /4 -B\\gf. 

Proof. Let us write g — (g + ,g~) T and gf :— g ± {-,j)- By Equation (|43|) we have 

\\4 2) gf n(2) = J2 (\\(dr m^ 1 + A(r))gff 

jez, 

+ \\(-d r -m J r- 1 +A(r))gJ\\ 2 ). 

Furthermore, dropping the term — d 2 , we get 

||(±9 r - m.r- 1 + A(r))gf\\ 2 > (gf | {{m 2 T m^ 2 + A{rf)gf) 

+ (gf | T d r A(r) - 2m J r" 1 A{r))gf) . 

Observe that (A2) implies that 



1 f r 

(75) -~ \ b(s)sds — o(l) , as r — > oo . 

r Jo 

This can be seen by splitting the integral above in the regions where b(s)s is inte- 
grable and the one where b decays in the L°°- norm. Hence, given q% £ (q±, 1) we 
find, using (|75|) . a constant R 2 > Ri such that, for all r > R2, 

B{r)>q 3 B , A{r)> qi B r/2, 

(76) 

\d r A(r)\<B, A{r)<Br/2. 

Therefore, for all r > R > R2, we get 

\\(±d r -m jr - l +A(r))gff 

> (df, (-^V4 + qlBlr 2 /A - 2m ] r- 1 A{r) - B)gf) , 

where we also use that (m 2 ± mf) > —1/4. 

Assume that rrij < 0. Since q\ > q 2 and A(r) > 0, for r > R2, we find that 

(77) \\{±d r - mj r- 1 + A{r))gf\\ 2 > (gf \ (q 2 B 2 r 2 /4 - r~ 2 /4 - B)gf) . 

Assume now that rrij > 0. Recall that A(r) < Br/2, for r > R 2 . Using that 
m, < r 2 {q\ - n 2 ql)B^/(4B) on supp# C £l qufiq2 we get 

W^dr-m^ + AirVgff 

> (gf I {q 2 B 2 r 2 /4 - r~ 2 /4 - m,B - B)gf) 

>(gf\^ 2 q 2 B 2 r 2 /A-r- 2 /A-B)gf). 

This finishes the proof. □ 



Proof of Lemma [3 Notice that 

e ±S K {2) e^ 5 = K (2) + Z ±p % Z ±p < := ±iSd r p e a 2 . 
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Thus, we have 

Q = Re((i^ 2) + Z-*)g | {K {2) + Z*)g) 
= \\K™g\\ 2 -5 2 \\d r p e g\\ 2 
Since \d r p € \ < \d r p\ < q 2 B r/2 we find 

(78) Q>\\K^g\\ 2 -{l/A)5 2 q 2 B 2 \\rg\\ 2 . 

Combining this with Lemma [TOl and that suppg C {(r,j) \ r > R} we obtain (recall 
that < 5 < p < 1) 

(79) Q > ([p 2 - S 2 )q 2 2 B 2 R 2 /A - R- 2 /4 §) \\g\\ 2 . 

This concludes the proof. □ 
Proof of Lemma [H We clearly have 

(80) Q<\{K^e s ^g,f R (E-W)i!}\ + \{K^e 5 ^g,a 2 (d r f R )^}\. 

(2) 

We analyze each of the above terms separately. Using that (Kq ' + W)^ = 
and noting that Wf R extends trivially to a bounded operator (for R > Ri large 
enough), we have, for any r\ £ JT7Qf(R 2 ; C 2 ), 

(K^r,, f R (E - W)4>) = ((E - W)f R K {2) V , $) 

= {K {2) f R (E - W)r), *) + ([(E - W)f R , K^]r,, ¥) 

= (t), (E - W) 2 f R $) + (r), [W,K^]f R ^) 

+ ( v , i <7 2 (d r f R )(w -E)fy. 

This identity extends to any r\ £ V(Kq ), in particular, we may choose r\ = e p 'g 
(see Lemma [7]). Thus, using Lemma [6j we find a constant C > 0, independent of 
R and e, such that 

\{K {2) e s r*g,f R {E-W)$)\ 

< \\g\\ \\e 5 "' [(E - W) 2 f R + [W, K (2) ]f R + i<J 2 (d r f R )(W - E)] || 

<C\\g\\\\e s ^<C\\g\\(\\g\\ + ^(1-^)11). 

We now treat the second term in (|80l) . We define the operators T and L acting, for 
any h £ H {2) and (r,j) £ R+ x Z, as 

{Th)(r,j) = e-^h(r,j), 

(Lh)(r,j) = (2a 1 a 2 (m j r- 1 + A(r))(d r f R Th)(r,j) . 

Clearly, since A(r) is bounded on the support of d r f R (for R > R\ large enough; 
see ([7610 L is an anti-symmetric bounded operator on 'H t ^ 2 \ With these definitions 
we have, using again the eigenvalue equation, that for any r\ £ J 7 UC^°(M 2 ; C 2 ) 

(K {2 \,a 2 {d r f R m 
= (K( 2) a2 (d r f R )r), *) + (t?, l S u PP d rfR (id 2 f R ^ - LT^ 1 ^)) 
= (V, \u PP drf R ^2(d r f R )(E - W)V +id 2 f R ^ - LT^ 1 *)) . 
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Note that T" 1 * e U {2) by LemmaE Next, we extend this identity to 77 € T>(K^ 2) ) 
and replace 77 by e Spc g. Using that e 5p£ l S upp a r f R is bounded uniformly in e > 0, 
we find e-independent constants C(R),C'(R) > such that 

\{KPe s o*g,a 2 (d r f R )4>)\ 

< C'(R)\\g\\ lle^ppa^lKHT- 1 ^! + ||l supp ^ /R iy*||) 

<C(R)\\g\\, 

where in the last inequality we use again Lemma [SJ Therefore, we obtain from (|80[) 
and the above bounds that 

Q< \\g\\{C\\g\\+C\\e s »(l-f R )\\+C(R)), 

which concludes the proof. □ 

Appendix A. Bounds for the Green function of D Ao 

Let 

0(x-x):= , t?(x,x):= — — (x 1 x 2 -x 2 x l ). 

Lemma 11. Let z 6 R\ct(Dao) an d ^ G*o(x, x', z), x, x' £ M 2 , be a representation 
of the Green kernel of (-Da ~ z) as 2 x 2-matrix. Then we have that 

(81) ||G (x,x';z)|| cw <e- e ( x - x Vx-x';z), 
for some function u>(-; z) : M 2 — > R + that satisfies 

(82) sup |x|e~ £|x| w(x; z) < 00, e > . 

x£R 2 

Proof. Recall that by Proposition [T] we have for E ^ that ±E £ o-(Dj± ) if and 
only if E 2 € a{dd*) \ {0} = a(d*d) \ {0}, where 

(83) d*d = (p - A ) 2 - B Q , dd* = (p - Ao) 2 + Bo . 
A simple computation using ([7]) yields, for any zfl \ a(Dji L0 ), 

(D^-zy^iD^+z^D^-z 2 )- 1 

(84) _ ( z{d*d~ z 2 )- 1 d*(dd*-z 2 )- x 



d(d*d-z 2 )- 1 zidd^z 2 )- 1 
It is well-known that the Green function of (p — Ao) 2 is given by 

(85) [(p - Ao) 2 - C] _1 (x,x') = (4 7 r)- 1 r(a)e i ' ?(x ' x ' ) e- e(x - x ' ) C/(a, 1, 20(x - x')) , 

where U is a confluent hypergeometric function and a = — l/2(£/i3o — 1) ^ — N, 
see for instance [5J Lemma 2.2]. 

Combining (|53")) . and (1551 we obtain that the Green kernel of Da is given 

by 

(86) Go(x,x';z) = e i ^ x ')- 9 ( x - x ') ( ^f^'; z \ ^(x^z) \ 

\ i2i 2 (x,x';z) S2 22 (x,x'; z) J 

where we define a± = —l/2((z 2 ± Bq)/Bq — 1) and 
ftii(x,x';z) := (47r)" 1 zr(a + )[/(a + ,l,26»(x-x')), 

ni 2 (x, x'; z) := (47r)- 1 B r(a_ + 1) U{ot- + 1, 2, 20(x - x')){i (a* - x[) + (z 2 - 4)} , 
ft 22 (x,x';z) := (47r)- 1 zr(a_)C/(a_,l,26»(x-x')). 
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Here we also used that -£U(al,t) = -ct-U(a- + l,2,t) (see [1, Eq. 13.4.22]. Since 
— a± $ No, the bounds (j8l"T) and ([82} follow now from the asymptotic formulas for 
U [JJ Eq. 13.5.2, Eq. 13.5.7, Eq. 13.5.9]. □ 

Appendix B. The family {H^(z)} ze c 

Throughout this section we assume that (A1)-(A4) are satisfied and use that no- 
tation introduced in Section[6l Our concern is the family of operators {H^ 1 ' (z)} 26 c 
defined a priori on the dense subspace UC§°(R 2 , C 2 ) of Ti^ as 

(87) H^(z) := K { 1} +v z , z£C. 

We first state a technical lemma. 

Lemma 12. Let T be a (complex-valued) multiplication operator on L 2 (R 2 ,C 2 ) 
with T £ Lf oc (]R 2 , C 2 ), p £ (2, oo] and lim„^.oo ||l{| x |> n }T'|| 0O = 0. Then, T is 
relative + compact. 

Proof. For n £ N write T = T x +T 2 where T x is supported inside the ball B n (0) C M 2 
and Ti on the complement of B n (0). Then T% is relative ■\/P 2 + l- compact [2"TI 
Theorem 4.1]. Moreover, 

||T(p 2 + l)- 1 / 2 -T 1 (p 2 + l)- 1 / 2 || < ||T a ||->0, 
as n -4 oo, from which follows the claim. □ 

Lemma 13. For any z € C £/ie operator v z (K^ + i) _1 is compact in HS 1 ^ . 

Proof. Let z e C and r > with r > |z|. Due to the inequality \v z \ < u T on 
M+ x T and the fact that u T £ L P (R + x T,rdrd0) (for some 2 < p < oo) we see 
that v z is well defined on the domain of Kq . Let " T = U*u T U . It suffices to show 
that U*v z U(Da_ + is compact in L 2 (R 2 ; C 2 ). This is, however, a consequence 
of Lemma and the discussion at the end of the proof of Lemma Q~] □ 

Lemma 14. {H^ 1 ' '(z)} ze c defined in (|87|l extends to an analytic family of type 
(A) with domain V(HW(z)) = V{K ( p). 

Proof. Due to Lemma [TBI we know that, for any z £ C, H^'(z) extends to a 
closed operator with T>{H^ 1 \z)) = T>(Kq). It is enough to show that, for any 
ip £ V(K { 1} ) the mapping C9z4 H^\z)tp £ is analytic. 

By the assumption (A5) we have, for any (r, 6) £ R + x T, that the power series 
«z(r,0)=£„ eNo vW(r,0>»with 



27ri jTj c|=s C" +1 

for some s > 0, has an infinite convergence radius. In addition, we clearly get from 
(|88f that |u( n )(r,0)| < u 2s (r,6)/s n for any (r, 0) el+xT. In particular, we find 
that 

||« (n Vll < ill wll , <pev(K^ 

Therefore, for any \z\ < s, 



s" 



neN 

This concludes the proof since s > can be chosen arbitrarily large. □ 
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